Three-state Potts model in combination with the rock-scissors-paper game 
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We study a three-state Potts model extended by allowing cyclic dominance between the states as it appears 
for the rock-scissors-paper game. Monte Carlo simulations are performed on a square lattice when varying 
the temperature and the strength of cyclic dominance. It is shown that the critical phase transition from the 
disordered state to the ordered one is destroyed by the cyclic dominance that yields a self-organizing pattern 
even at low temperatures. The differences and similarities are discussed between the present model and the 
half-filled, driven lattice gases with repulsive interaction. 
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The ordering phenomena and the related phase transitions 
are already well understood in the equilibrium systems |1] 
while the theoretical understanding of the non-equilibrium 
phase transitions is still at its beginning (3 1^ ■ Many rele- 
vant and general features of these transitions in the equilib- 
rium systems can be studied by the Potts models f7|. Now we 
introduce an extended version of the three-state Potts model 
to investigate the effect of cyclic dominance between the 
states. This model can be considered as a combination of 
the traditional Potts model and a spatial rock-scissors-paper 
game (sometimes called as three-state cyclic predator-prey or 
Lotka-Volterra models). The strength of cyclic dominance 
will be characterized by a single parameter (e) in such a way 
that for e = the system becomes equivalent to the equilib- 
rium Potts model exhibiting a well-known critical transition. 
The consideration of this model was strongly motivated by the 
work of Katz et al. 1 5, 6] who introduced the concept of driven 
lattice gases to study the effect of an external electric field on 
the ordering process. Since the appearance of their pioneering 
works many general features of these systems have already 
been explored (for a review see the Refs. |2, 7]). 

The present model provides a continuous transition be- 
tween an equilibrium system and a spatial evolutionary game 
where the time-reversal symmetry (detailed balance) is broken 
at the elementary (microscopic) steps. It will be shown that in 
the presence of cyclic dominance the long-range order cannot 
be observed even for low temperatures, and thereby the crit- 
ical transition is also suppressed. A similar phenomenon has 
already been observed for the two-dimensional driven lattice 
gas with repulsive interactions where the formation of long- 
range order is prevented by an interfacial instability due to the 
enhanced particle transport along the boundaries separating 
the "chessboard" and "anti-chessboard" ordered phases |^|^. 

Our analysis is focused on a two-dimensional system where 
each site a; of a square lattice is characterized by a three-state 
site variable, namely, s{x) — sq, si, and S2- For later conve- 
nience these states (strategies, species, etc.) will be denoted 
by the basis vectors of a three-dimensional space, i.e. 



The time evolution of the system is governed by random se- 
quential updates. More precisely, the transition probability 
from a state s{x) (at site x) to a randomly chosen state s'{x) 
is given as 



W\s{x) 



1 + cxp(-(5[/(x)/T) 



(2) 



where 5U [x) is the difference of payoffs between the final and 
initial states, and T is the temperature characterising the effect 
of the noise. The payoff at site x depends on s{x) as well as 
on the neighboring states [s(y)] as given by the following sum 
of matrix products: 
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where the summation runs over the nearest-neighbors of the 
site X, s+(x) denotes the transpose of s{x), and the payoff 
matrix A is defined as 



(4) 



In the limit £ — > this model can be considered as a (kinetic) 
three-state ferromagnetic Potts model with a Glauber dy- 
namics 1 10]. Evidently, in this case the total energy is defined 
as H ~ — U{x)/2 and the microscopic processes satisfy 
the detailed balance in equilibrium. Consequently, the system 
tends towards a stationary state whose statistical features are 
described by the Gibbs ensemble. When decreasing the tem- 
perature the Potts model undergoes an ordering process from 
the disordered state to one of the three equivalent homoge- 
neous (ordered) states. The correspondingcritical transition 
represents a well-known universality class 

For £ > the off-diagonal components of the payoff matrix 
A are asymmetrical therefore the total payoff (or the above de- 
fined H) is not affected by the value of e for any states. At the 
same time, the value of e influences the probability of strat- 
egy changes because VF[s(a;) s'{x)] depends definitely on 
the variation of individual payoff [SU (x)] and the above evo- 
lutionary rule manifests a way how the (selfish) individuals 



2 



wish to maximize their own payoff without any concern about 
the neighbors' performance. As a resuh, cycHc invasions oc- 
cur along the boundaries separating homogeneous domains; 
domains of state sq are invaded by si invaded by S2 invaded 
by Sq. These cycUc invasions are capable to maintain a self- 
organizing pattern with rotating spiral arms whose "velocity" 
is controlled by e fll[Tl[ll. 

For most of the spatial evolutionary games the choice of the 
dynamical rules (or M^[s(a:) s'{x)]) is based on a learn- 
ing mechanism or strategy adoption modelling the Darwinian 
selections 1 14. ITsl IT6ll . In these models different ways are 
suggested for the players to adopt the strategy of their more 
successful neighbors. The common feature of these strat- 
egy adoptions is that the new state will be equivalent to one 
of the neighboring strategies. Consequently, this mechanism 
prohibits the variation inside the homogeneous domains and 
makes the extinction process to be similar to those defined 
by the contact process (or directed percolation) f^, '3']. In the 
present model, however, the above "Glauber dynamics" al- 
lows the players to choose all the possible strategies therefore 
the time variation is not restricted to the interfaces separat- 
ing the homogeneous domains. In the context of evolutionary 
game theory the above evolutionary rule describes a different 
behavior Namely, here the players know all their possibilities 
and their choices depend on the increase of income what they 
are able to evaluate in the knowledge of neighboring strate- 
gies. 

In the present work we study the effect of the "cyclic dom- 
inance" on the phase transition. For this purpose systematic 
Monte Carlo (MC) simulations are performed on a square lat- 
tice under periodic boundary conditions varying the tempera- 
ture T and strength e of the cyclic dominance for different lin- 
ear sizes L. Each simulation is started from a random initial 
state and after a suitable thermalization time we have recorded 
the concentration of states (pQ, pi, and P2) for each Monte 
Carlo steps. We have also made simulations starting from or- 
dered homogeneous phases to check the stability of the sta- 
tionary state. To investigate the ordering process, we have 
determined the average value of the order parameter from the 
values of concentration data 



tion that deviates monotonously if we decrease the system size 
1,17.. 18 1. Significantly different finite size effects are observed 
when investigating the present model for e > 0. 
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[3max(po,Pi,P2) - 1]) 
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where (• • •) refers to averaging over a sampling time var- 
ied from 10^ to 10'' Monte Carlo steps per sites (MCS). In 
the disordered phases po = Pi = P2 = 1/3 (due to the 
cyclic symmetry) and m = in the thermodynamic limit 
(L 00). For e = and below the critical tempera- 
ture (T < Tc = 0.995(1)) the system evolves into one of 
the long-range ordered (symmetry breaking) stationary states 
(e.g., (po) = (1 + 2to)/3 and (pi) (^2) = (1 - m)/3 and 
the remaining two equivalent states are given by the cyclic 
permutation of indices) if the linear size is sufficiently large. 
In the thermodynamic limit the order parameter m decreases 
monotonously from 1 to as the temperature is increased from 
to Tc and the vanishing of m follows a power law behavior 
if Tc is approached from below |4]. For finite sizes, however, 
the MC simulations exhibit a smoothed order parameter func- 




FIG. 1 : Monte Carlo data for the order parameter vs. temperature at 
£ = 0.1 for different system sizes as indicated The solid lines are 
guides to the eye. 

Figure ^illustrates how the order parameter m varies with 
temperature T for different linear sizes if e — 0.1. Apparently 
the MC data refers to an ordering process for small sizes (L = 
50 and 100) bearing a resemblance to MC data obtained for 
e = 0. On the contrary, for example when L = 500, the 
MC data do not indicate the appearance of long-range order. 
Instead of it a self-organizing, three-color domain structure 
can be observed when visualizing the time-dependence of the 
spatial distribution (for a snapshot, see Fig.|3- 




FIG. 2: Snapshot on a typical domain structure appearing for T — 
0.64 and e — 0.1. The three-edge vortices (antivortices) rotate in 
clockwise (anticlockwise) direction with spiral arms because the av- 
erage velocity of the invasion fronts (white invades black invades 
gray invades white) is hardly affected by their curvature. 

On this snapshot one can identify all the three ordered 
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phases forming domains with a characteristic linear size I. For 
T < Tc and e = the growth of these domains (7~ \/t) is 
driven by the interfacial (Potts) energy IT^ Eol I21I1 . Here, 
however, this domain growth is prevented by two processes 
emerging for e > 0. The first process is related to the appear- 
ance of rotating spiral arms for the three-edge vortices where 
the three types of domains (or domain boundaries) meet. On 
these maps we can distinguish vortices and antivortices rotat- 
ing in opposite directions. Some topological and geometrical 
features of such spatio-temporal patterns were already inves- 
tigated in previous papers 1 1 1, 12, 13, 22]. It is found that the 
spirals become well-marked if a "surface tension" is switched 
on and then the corresponding patterns cannot be character- 
ized by a single parameter (e.g., typical domain size or corre- 
lation length) II22I1 . 

In the present model there exists a second process causing 
the appearance of growing domains via a nucleation mecha- 
nism inside the large "homogeneous" territories. Due to this 
process an "ordered state" prevailed by sq will be transformed 
into another one prevailed by si as indicated by simulations 
for small sizes. In these cases the three ordered states fol- 
low each other cyclically and the above method yields a suffi- 
ciently large value for m (see Fig.Q. Both the duration time 
and the probability of these transitions, that are initiated by 
a nucleation mechanism due to the thermal fluctuations, in- 
crease with the system size. This is the reason why the values 
of m are higher for L = 50 than those for L = 100 in the 

Fig.m 

For sufficiently large system sizes (L >> I) both men- 
tioned mechanisms work simultaneously and result in a self- 
organizing pattern where the three states are present with the 
same concentration and m ~ 0). Henceforth the quantitative 
investigations will be focused on the large systems {L > 500) 
and to the region of temperature (T > 0.6Tc) where the spa- 
tial patterns are isotropic. 

Now we study the variation of correlation length ^ derived 
from the asymptotic behavior (exponential vanishing) of the 
two-site correlation function f§\. For this purpose a series of 
MC simulations is performed by varying the temperature for 
e — 0.1. The inset in Fig.|3lillustrates the absence of diver- 
gence in ^ as expected. When decreasing the temperature the 
correlation length increases monotonously until a maximum 
value (^ ~ 11.5(5) measured in lattice unit). Below the peak 
at T ~ 0.77 the visualization of the distribution of species 
shows a self-organizing pattern (see Fig.O and here the value 
of ^ decreases very slowly with T. 

The e-dependence of the correlation length is also deter- 
mined for a fixed temperature and the results are illustrated 
in a log-log plot (see Fig. |3}. These MC data are consistent 
with a prediction ^ 1/e for small e. In the typical size of 
domains similar divergence was found previously for a model 
where the nucleation mechanism was blocked 12?]. This ob- 
servation refers to a minor role of nucleation mechanism in the 
maintenance of the self-organizing patterns at sufficiently low 
temperature {T < Tc). For higher temperature the nucleation 
mechanism plays a crucial role by preventing the formation of 
monodomain state even for e = 0. Thus it is conjectured that 
this process results in a different behavior of ^ at the critical 
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FIG. 3: Variation of correlation length (with error bars) as a function 
of e for fixed temperature (T — 0.75) The straight line indicates a 
power law divergence with a slope of — 1. Inset shows the correlation 
length vs. temperature for e — 0.1. 



temperature. Furthermore, here it is worth mentioning that in 
the above mentioned driven lattice gas model the transversal 
correlation length was also proportional to the inverse of the 
strength of driving field . 

In such systems the Potts energy measures the concentra- 
tion of domain walls and it gives an additional information 
about spatial distributions. From the average Potts energy 
as a function of temperature one can derive a specific heat 
(c — d{H)/dT) that exhibits a A-divergence at the critical 
temperature in the equilibrium limit (e = 0). Figure 0] il- 
lustrates how the A-divergence is smoothed out if the cyclic 
dominance is switched on. When choosing larger and larger 
e the maximum value of specific heat decreases meanwhile 
the peak position moves towards the lower and lower temper- 
atures. The peaks are so shallow in the "driven" cases that a 
logarithmic scale was necessary to present them in the same 
figure. The appearance of the this peak in the specific heat can 
be interpreted as a sign of the short range ordering. 

A similar phenomenon was observed for the driven lattice 
gases with repulsive interaction when increasing the external 
electric field |8, 9]. Although the observed patterns and mi- 
croscopic mechanisms are very different, in both cases inter- 
facial effects prevent the formation of long-range order in the 
presence of driving force. In these cases the interfaces belong 
to the stationary states and their geometrical characterization 
requkes additional parameters. This general feature can oc- 
cur in some other non-equilibrium systems (e.g., in ecological 
models) where an external force induces some extra activity 
along the interfaces separating the "ordered domains". 

In summary, a three-state dynamical lattice model is intro- 
duced by combining the Potts model and the rock-scissors- 
paper game to study the effect of cyclic dominance on the or- 
dering process. Due to the cyclic dominance the time-reversal 
symmetry is broken at the elementary steps and thereby the 
behavior of this this model cannot be described by the meth- 
ods of equilibrium statistical physics. Our numerical analyses 
have justified that both the long-range (symmetry breaking) 
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FIG. 4: Specific heat as a function of temperatiire for tiiree different 
values of e as indicated. These MC data are obtained for such a large 
Unear size (L = 1000) where the size effects are already negligible. 



ordering process and the corresponding critical transition are 
suppressed in the presence of cyclic dominance (e > 0). Ac- 
cording to the simulations the three equivalent ordered phases 
coexist by forming a self-organizing domain structure even 
at low temperatures and sufficiently weak cyclic dominance. 
The equilibrium state is approached via the divergence of the 
typical domain size when the strength of cyclic dominance 
goes to zero. 
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